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TECHNICAL FIELD OF THE INVENTION 

This invention relates to the design of optimal multivariable linear controllers and PID 
controllers. 

BACKGROUND OF THE INVNTION 

A PID (proportional-integral-derivative) controller is a device, usually implemented in a 
computer, that is used to control a process variable (e.g., temperature, pressure, etc) of an 
industrial process. For a multiple-input and multiple-output (MIMO) process, the process 
variable y(t) and controller output u(t) are n-dimensional vector variables of the time t The 
controller receives the sample values of y(t) and calculates its output u(t) according to a PID 
equation. The controller sends its output u(t) to the process so that the error e(t)=r(t)-y(t) 
approaches zero as time t increases, where r(t) is the set point (also known as the 
reference signal or the command signal, etc). 

There are many types of PID controllers, depending on the use of different types of PID 
equations. In discrete time form, all these PID controllers can be viewed as special cases of 
the general linear controller with the linear control equation 

Du k =Er k -Cy k 

where y k =y(tk) is the process variable at time t k =to+kT 8t to is the initial time, T s >0 is the 
constant sampling period, k=0, 1, 2, ... is a non-negative integer called discrete time 
variable, u k =u(t k ) is the controller output at time t k , and u k can be subject to upper limit 
and/or lower limit constraints placed on one or more of its components, r k =r(t k ) is the set 
point at time t k , D, E and C are n by n matrix polynomials in the unit backward shifting 
operator z' 1 such that for any discrete time signal x k , z 1 x k =x M , and one or more of the said 
D, E and C contain tuning parameters that are to be determined, 
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For example, if D=(1- z" 1 )/T s *l, where I is the identity matrix of order n, and 
E=C=Ki+KpD+K d D 2 , where Kj, K p and K d are user specified constant coefficient matrices 
called the integral gain, the proportional gain and the derivative gain, then the linear 
controller becomes the following "type A" PID controller 

Du k =Kie k +K p De k +KdD 2 e k 

where e k =r k -y k is the error. The three gains K it K p and Kd are the tuning parameters. 

If the aforesaid E is changed to Kj+K P D or Kj, then the linear controller becomes the 
following "type B" and "type C" PID controllers, respectively: 

Du k =Kie k +K p De k -KdD 2 y k 
Du k =Kie k -K p Dy k -KdD 2 y k 

Once its structure is selected, the performance of a linear controller or a PID controller 
depends mainly on the choice of its tuning parameters. 

How to properly choose the tuning parameters for a PID controller is a problem that has 
attracted a lot of studies ever since PID controllers became widely used in industry in the 
early 1940s. The Ziegler-Nichols tuning methods developed by Ziegler and Nichols in 1942 
and 1943 are still widely used in industry. Other model based methods choose the tuning 
parameters by minimizing some well-known control performance index such as the integral 
absolute errors (IAE), the integral squared errors (ISE) t etc. However, practice shows that 
all these methods often lead to oscillatory control results. 

DETAILED DESCRIPTION OF THE INVENTION 

Suppose the open-loop transfer function of the process from u k to y k is A" 1 B, where A and B 
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are known n by n matrix polynomials in the unit backward shifting operator z~ 1 : 

A=l-A 1 z 1 -A 2 z" 2 -... -A p z p 
B= BiZ' 1 + B 2 z 2 + ... + B q z" q 

where I is the identity matrix of order n, p and q are two known positive numbers, and A 1t 
A 2 , .... A p , Bi, B 2l and B q are constant coefficient matrices of order n. A and B can be 
obtained by use of the commercially available software product "System Identification 
Toolbox" from The MathWorks Inc. Below is another way for finding A and B. 

Select a positive integer v that is not less than n(p+q) and let k be not less than v + max{p, 
q}, where max{p, q} denotes the larger one of the p and q. Let h k be a vector defined by 

h k T= [yk-i T y k -2 T .-. yk- P T u k .i T u k _ 2 T ... u k . q T 1] 
where for any matrix X, X T denotes the transpose of X, and let 

Y=[y k y k -i ... y k -v] T 
H=[h k h k .i ... h k _ v ] T 
W=diag(1, w, w 2 , .... w v ) 

where w is a positive real number not larger than one, and W is a diagonal matrix with 1 , w, 
w 2 , and w v on its diagonal. The number k should be large enough and the controller 
output signal u k should be chosen so as to ensure that the matrix H is of full column rank. 
Then the matrices A!, A 2| .... A p , Bi, B 2l .... and B q , and the vector d can be obtained from 

[Ai A 2 ...A P Bi B 2 ... B q d] = ((WH) + (WY)) T 
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where for any matrix X, X* denotes the Moore-Penrose pseudo-inverse of X, which is a 
standard function under the name "pinv" in the software product MATLAB from The 
MathWorks Inc. (MATLAB is a registered trademark of The MathWorks Inc.). Thus the two 
polynomials A and B and a vector d are obtained. When A t B, and d are so obtained, the 
process can be approximately represented by the following linear model: 

Ay k =Bu k +d 

And the open-loop transfer function of the process is then A" 1 B. Please note that this 
invention does not deal with the problem of how to get the polynomials A and B. From now 
on it is always assumed that, one way or another, the two polynomials A and B are known. 

If the process is controlled by the following linear controller as described in the previous 
section: 

Du k =Er k -Cy k 

where D, E, and C are n by n polynomials in the backward shifting operator z \ and one or 
more of D, E and C contain tuning parameters that are to be specified, then the closed-loop 
transfer function from the set point r k to the process variable y k is 

Q=(A+BD' 1 C)- 1 BD" 1 E 

From Q the characteristic polynomial in z, denoted by b(z), can be found. All roots of the 
characteristic equation b(z)=0, denoted by z 1t z 2 , z 3 , Zj, where j is a positive number, 
form all the poles of Q. 

Let max {(z^, |z 2 |, |Zj|} denote the maximum of all the absolute values |zi|, |z 2 |, and 



